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ABSTRACT 


The theoretical formulas for the influence of the change of moment of inertia and mass-loss on the variation of the pulse period in 
the magnetic dipole model are given. The numerical solutions of slow down of period due to mass-loss are calculated by using the 
given formulas for three pulsars: PSRO355+54, PSR1930+22 and PSR1800-21. The numerical results are listed in Table 1. Discussion 


and conclusion are drawn. 
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INTRODUCTION 


Some authors studied the variation of Pulse Period arisen from the change of moment of inertia. However, they always use the 
method of the angular momentum conservation (IO= const). The angular momentum is not conservative, if when we consider 
mass-loss of a star. Hence the change of angular momentum arisen from variation of mass-loss is not suitable to the research by 
using the angular momentum conservation. In the formula of the magnetic dipole model of pulsars the moment of inertia is not 
variable but when we consider the mass-loss of pulsar, the moment of inertia may be changed. Hence it is necessary to give the 
formula for the magnetic dipole model suited to the variation of moment of inertia. This is an important work in this paper. 

There are some authors for researching the influence of mass-loss of a star on the rotation .However, there is a few author for 
researching the influence of mass-loss of pulsar on the rotation. Esposito and Harrison (1975) researched the influence of mass-loss 
on the orbital period for the Taylor binary pulsar, but they did not research the influence of mass-loss through the change of 
moment of inertia on variation of the period. This paper studies some works in this aspect. 
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THE FORMULAE FOR THE INFLUENCE OF CHANGE OF MOMENT OF INERTIA ON THE 
VARIATION OF THE PERIOD OF PULSAR 


The pulsar radiating power W is transformed from the rotational energy at a rate ae P 
t 


dE dE 
W =-— or W+—=0 (1) 
dt dt 
According to the theory of magnetic model (Ostriker &Gunn,1969,and Shapiro & Teukolsky 1983) 
dE Z 52 20 
— =-—(M sina)’ Q* =-—— a , Q=— (2) 
dt 3c 3c OP Pr 
Substituting the equation (2) into the equation (1), we can obtain the radiating power W 
322° 
re 6) 
oe 


Where P is the period of pulsar, and wu is the projection of the magnetic dipole moment M on the direction perpendicular to the 


rotational axis. Q denotes magnetic inclination 

. . = ? 
uw” =(M sina)’ =(M,sina)’e at's go, (4) 
& is the coefficient of magnetic decay, 


The energy carried away by radiation from the rotational energy of pulsar can be written 
1 2 

B= 10? = 22°17, (5) 
2 P 


Here I is moment of inertia. If we consider the variation of moment of inertia with time, then 


1 dl 21(t) dP 
2p a, (6) 


Substituting the formula (3) and (6) into the formula (1), we obtain the Bernoulli equation for n=1 


dP 11d pe 8°" (1) 5 


(7) 


dt 2 1 dt 3c? I(t) 
We can transform Bernoulli equation into the first order linear differential equation. Both sides of the equation (7): are multiplied by 
2P. ie. 
dP I 162° yu? 
2P Pp? = lacy (7) 
dt I 3c°I(t) 
The equation (7) may be written as the form of the first order linear differential equation 
dP* I 
—~—P* =Q(t), (8) 
dt I 


I 
We define N = as 


According to the first linear differential equation (6), N is the function of time f or it is a constant value. In this paper N is a 
constant value as shown in the expressions (22). Integrating (9), one yields 

-Nt 
T=I,he”, (10) 
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Q(t) = 1677 u(t)? = 1627 uo eet 


3c71(t) 3c71, 
The equation (8) may be written as 
ip? 
2 + NP* =O(t). (11) 
t 


Integrating the equation (11), one yields 


PO)? =e I] J OG! #2: 


Substituting (10) into the above ide expression, we obtain 
167" Uo 
When we take t=0, P?(t) = on “C= PO)’ i.e 
P(t)? =e “[P()* + Loatte" } ‘eG 2N) ay] 
3c°I, °° 


Integrating the above expression, we obtain 


167° I 
P 2 _ -kt P 22: 0 ~(E-2N)t ash. 
(t)" =e [Po oa E*IN) a? ~) le } (12) 


When we only consider magnetic decay and do not consider the variation of moment of inertia, i. e 


N =0,,. € <0, 


—}). (13) 
oe Hy 


When we only consider the variation of moment of inertia and do not consider the magnetic decay, i. e. N #0, S = 0, then 


8r- I, : 
Se ay), (14) 


P(t)? =e “[P(0)* = 
() PO) = soy 


dl — 
In the equation (7) I = Io, a 0, P=P(0), P = P(O), as t =O we obtain 
t 


ee = PO) PO) 05) 
Substituting (15) into equations (13) and (14), the equations (13) and (14) can be written as 
P(t)? = P(O)* —2P(0)P(O\(e* —D/E (16) 
P(t)? =e {P(0)? = PO) P(O)[e2™ —1]}/N (17) 


Hence, we can estimate the variable rate of pulse period per century 


OP =|P(t)- P(t, )|(5/,_). (18) 


cent 
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Here P(0) is the initial value as t = 0. 


THE INFLUENCE OF MASS-LOSS OF PULSAR ON THE VARIATION OF ITS PERIOD 


dP 
Esposito and Harrison (1975) wrote the luminosity of pulsar B = 47mR : re, /5P* which equal to the rate of mass-loss- 
t 
Cc” —— . This is obtained from moment of inertia as a constant. In the present paper moment of inertia, I, is variable with time. We 


2 
use the formula: J =—mR - , where R is radius of pulsar which is assumed to be constant. 


Because mass is variable with time .Hence 

T 2 p) Ion 
cs R? ia ak. (19) 
dt 5 dt 5 I om 


dm 
Substituting these into the expression (5), and letting it equals — room ae , and then, we have 


e dm 4 2(R)dm __ 8 > mR° dP 


‘/s 
dt 5 P° at 5 P? dt 
ldm 8a°R’ dP, 42°R’., 82°R°’ dP. 42°R* 
Es =— ; (c* — S) => (l=) 20) 
m dt 5P- dt 5P 5c°P° dt 5c°P 


The second term of right hand side may be neglected because va is very small. According to the expression (9) and (19)-(20) we 
c 


can write 
I 1dm_82’R’ aP 


N= = =p ps 
I mdt 5c P° dt 


(21)= 


and P(t) 3 is variable with time.. 


Based on the formula (9), N is a constant, but in the right hand side of the above formula 


dP(t) 
dt 


dP(t) — dP(O 
a) 7 of) = const and P(t)’ = PO)’ = const. We may use Taylor series to expand 


In the following we may prove 


dt dt 
dP(t) and P(t) as follows 
t 
— = P(t) = P(O) + P(O)(t ty) + eveseeees 


P(t)? = PO)? +3P(O)? P(O)(t —t) + .oeeoeeees 

For a lot of pulsars P(O) ~ 0.1, P(O) ~ 107'°, P(O) ~ 10 we take t —t, = 100 year(Century) ~ 10°(s) 
P(O)(t—t,) ~ 10°" 
3P(o)’ P(O)(t-t,) ~ 10° 


=const and 


The terms P(O)(t = to) and 3P(0) 2 P(0) so small that may be neglected. Hence a = ae) 


dt 


P(t)? = P(O)* = constant. Therefore, in the formula (16) N is a constant. The exponential form (8) is hold well 
Pp 
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we I 1ldm_8 2°R* dP(O0 
I mdt 5c’P(0) dt 


= 2.5266x10° Ey) = Constant. (22) 
P(O) 


NUMERICAL RESULTS 


We only consider the case of the variation of moment of inertia due to mass-loss and do not consider magnetic decay. We use the 
formulas (17)- (18) and (22) to calculate the variation of pulse period due to mass-loss for PSRO355+54, PSR1930+22 and PSR1800-21. 


dP 
For these pulsars we assume that R = 12 x10°km. The data of P and ae of PSR0355+54 are adopted from Table of pulsar 
t 


dP 
parameters given by Manchester & Taylor (1977). The data of P and —— PSR1930+22 and PSR1800-21 are adopted from given 
t 


by Arzoumanian et al (1994). Substituting the above data into the formula (17) and (18), we obtain the numerical results for the rates 


of spin down of pulse periods due to mass-loss for three pulsars are listed in Table 1. 


Table 1 Numerical results of three pulsars 


P(0)/dt 
po) @ Old iat Didi. P(t) oP 
Pulsars 1075 N=- =— (l0”/s) , ‘ 
@ dos’) ce aa (Gf) (fF) 
S cent cent 
PSRO355+54 =—-0.1563 4.39 .90x102 0.1576 0.0013 
PSR1930+22 0.1444 63 5.28x1071 0.1463 0.0019 
PSR1800-21 0.1336 134.229 1.4225 0.1376 0.0040 


DISCUSSION 


(1) In the quadrupole elastic energy model of neutron stars the total energy and moment of inertia connects with oblateness & : 
1 2 2 2 
E=E, aoe +Ae°+Ble-é), 1 =1,d +6). 
But in the magnetic dipole model we may not consider oblateness or € = 0 . We consider pulsars as spherical stars.. 
2 -& . é : : ' fot ssa : 
(2) The formula LZ a My e * cited in this paper is suitable to the magnetic inclination a, as a constant. Because 
. 2 -& 3: 203 = De 2 20: 20: : 
uw =M’sin*a=M,e"sin’a=(M, sin’ae“=",e" My =M, sin*>a=M, sin’ aie, 
a =A, =const, 
Hence the research of this paper includes the magnetic inclination as a constant. However the magnetic inclination of pulsars 


varies very small in a century. We can assume that the magnetic inclination a is a constant in a century (3). We may discuss mass- 


less as compared with magnetic decay. We use the formula (16) to calculate the variable rate of pulse period due to magnetic decay. 


Substitution of these data € = 1.3 x 10°/ yror & =3.54x 10°*/s given by Qu et al (1976) into the formula (16), we obtain, 


For PSR0355+54: (OP) = 0.00029( y cent ), 


mag 


For PSR1930+22: (OP) ne 7 0.000327), 


For PSR1800-21: (OP) mag = 0.00034( vA vr p- 


From these numerical results as compared with the results due to mass loss in Table 1 we can see that the values for magnetic 
decay are nearly the same order with that of mass-loss for three pulsars. Hence both results cannot be ignored. 
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CONCLUSIONS 


(1) The change of moment inertia of pulsars can influences spin down of pulse period in the magnetic dipole model, and spin down 


connect with the exponential formulation. 


(2) Pulsar loses its mass with the formulation of the magnetic dipole radiation. It can influences spin down of pulse period through 
the change of moment of inertia. The rate of spin down of pulse period is the order 10° cae Py for three pulsars. The rate of 


slow down due to mass-loss can be observed by using the recent astronomical instruments through a long time. 
(3) It can be seen from Table 1 that the shorter the pulse period, the longer the rate of pulse period per century, such as PSR1800-21. 
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